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In the Hamiltonian approach to Coulomb gauge Yang-Mills theory, the functional Schrodinger 
equation is solved variationally resulting in a set of coupled Dyson-Schwinger equations. These 
equations are solved self-consistently in the subcritical regime defined by infrared finite form factors. 
It is shown that the Dyson-Schwinger equation for the Coulomb form factor fails to have a solution 
in the critical regime where all form factors have infrared divergent power laws. 
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I. INTRODUCTION 



Low energy gluon modes are expected to be responsible 
for the distinctive features of QCD such as confinement or 
chiral symmetry breaking. Current evidence for gluonic 
excitations is sparse 0, US 3 > however new experimen- 
tal efforts at JLab, PANDA, and BES, through studies 
of hybrid and glucball spectra are expected to shed more 
light into the nature of gluonic excitations. It is thus 
highly desirable to develop a theoretical framework for 
studies of the Yang-Mills sector that is both: rooted in 
QCD and practical. Coulomb gauge q uantization offers 
such a framework @, @, 0, 1, H, M, [Li El El ■ Elimina- 
tion of the longitudinal component of the gauge field by 
imposing the Coulomb gauge condition, V • A a (x) where 
a denotes color (a = 1 • • ■ Nq — 1), leaves only two physi- 
cal transverse components as independent degrees of free- 
dom. In the functional integral approach, the Aq compo- 
nent becomes constrained and defines the instantaneous 
potential of color charge, while in the canonical quanti- 
zation approach Wcyl gauge, Aq = 0, is used and the 
static potential emerges from the resolution of the Gauss 
law. In the Coulomb gauge the resulting many-body 
system of transverse gluons interacting via long-range 
Coulomb exchange forces can be studied using standard 
techniques. After ultraviolet (UV) regularization it is 
possible to introduce an ansatz for the ground state vac- 
uum wave functional in terms of the dynamical gluon 
variables and optimize it by varying the energy density. 
In particular with a gaussian ansatz a canonical trans- 
formation exists which transforms the gauge fields to the 
particle representation. These quasi-particles satisfy a 
dispersion relation with a mass gap that originates from 
non-pcrturbative self-interactions mediated by the long 
range Coulomb potential. The expectation value of the 
Coulomb potential is computed self-consistently in the 
same vacuum state, and these self-interactions lead to a 



strong enhancement of the potential for large separation 
between localized color charges. For a true linear poten- 
tial the single quasi-gluon mass becomes infinite which is 
consistent with what is expected for a confining phase: 
formation of color states requires infinite energy. In a 
color singlet combination, however, residual interaction 
between constituents screens the long range interactions 
and leads to a finite mass for bound states. A priori, 
however it is not guaranteed that a particular ansatz 
for the vacuum wave functional would result not only 
in qualitative but also in quantitative description of all 
features of confinements. The gaussian ansatz represents 
fluctuations around topologically trivial distributions of 
the gauge field and since topologically nontrivial config- 
urations can also minimize the non-abelian Yang-Mills 
action, a more complicated vacuum wave functional may 
be needed to, e.g., reproduce the area law behavior of 
the spacial Wilson loop [3, EH- Also, the vacuum ex- 
pectation value of the Coulomb potential is actually not 
the same as the energy of the state with static sources. 
External sources polarize vacuum and lead to string for- 
mation for large separations. Thus one should not com- 
pare the expectation value of the Coulomb energy to the 
"Coulomb plus linear" potential from lattice gauge sim- 
ulations of temporal Wilson loops [H, E3] ■ Nevertheless, 
the variational Coulomb gauge approach does reproduce 
the qualitative features of confinement and studying dis- 
crepancies between this approach and lattice results, can 
improve our understanding of the underlying many-body 
problem. 



In a scries of papers we studied solutions of the varia- 
tional problem for the vacuum and single quasi- par ticle 
properties under various approximations [n], [Til. Il2l. El, 
EsLLSl]- The goal of this paper is to review these results, 
clarify the differences of the various approaches and re- 
solve some open problems. Ref. [TcJ EH and [H, EH 
use different ansatze for the wave functions and differ 
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in the extend to which the curvature in the space of 
gauge orbits introduced by Coulomb gauge fixing was 
included. While it was shown that the different ansatze 
for the variational wave functionals is irrelevant to the or- 
der of approximation [l3j , in particular for the infrared 
behaviour, the inclusion of the curvature of the space 
of gauge orbits is crucial for obtaining the correct in- 
frared behaviour. The subject of the present paper is 
the following: We reconsider the renormalization of the 
Dyson-Schwinger (DS) equations resulting in the varia- 
tional approach. In particular, we show that to the order 
considered all UV-divergencies can be removed by adding 
appropriate counter terms to the Hamiltonian. While in 
previous calculations the horizon condition (i.e. an in- 
frared diverging ghost form factor) was either assumed 
or at least seem to follow from numerical solutions in the 
present paper we study the coupled DS equations in the 
subcritical regime abandoning the horizon condition. In 
particular, we will use the infrared analysis of Ref. (l9| 
to investigate the criticality of the various DS equations, 
i.e. the disappearance of the solution of the DS equations 
as the infrared exponents exceed certain critical values. 
We will find that in the limit of critical coupling where 
the ghost form factor becomes infrared singular, the DS 
equation for the Coulomb form factor ceases to have a 
solution. This explains why no consistent solution for 
the Coulomb form factor was found when the horizon 
condition is implemented. 

The organization of the paper is as follows: In the fol- 
lowing Section we summarize the approximations used 
in solving Dyson-Schwinger equations of the variational 
Coulomb gauge problem from earlier analyses. In Sec- 
tion IIIII we derive the DS equations and discuss various 
renormalization schemes. In Section IIVI we discuss the 
infrared (IR) limit of the equations. Numerical results 
and discussion is presented in Section [V] Conclusions 
and outlook are given in Section IVI1 



II. VARIATIONAL COULOMB QCD 

The set of Dyson-Schwinger (DS) equations for correla- 
tion functions describing gluon, ghost and the Coulomb 
propagators are considered within a variational ap- 
proximation for the ground state vacuum wave func- 
tional. The Yang-Mills Coulomb gauge Hamiltonian 
H = H(A, II) is a function of the transverse gluon field, 
A a (x) and the conjugated momenta, Il a (x) = —E a (x), 



[A a (x),U a (y)} 



iS ab 5 T (x 



y) 



(i) 



with St = (1 — W/V 2 )5. In the Schrodinger repre- 
sentation, matrix elements of an operator C(A, II) are 
given by 



(*'|C>|*) = / VAJ{A)y'(A)0 



A l 5A 



*{A). (2) 



In [I(J, [IH the variational wave functional ^(A) = 
(A\Q) = ^o^O; A) m form of a gaussian ansatz was used, 

*o(^n,^) = cx p(-^ J dxdyA a {x)ujo{\x - y\)A a (y) 
expf-J / 7^A a (k)tJ (k)A a (k) 



2 J (27r) 3 ' 
cxp (~ J AujoA^j 



(3) 



where A a (k) is the Fourier transform of A a (x). The 
functional integration measure contains the Faddeev- 
Popov (FP) determinant, J = expTr ln(— Dd) where 
D(x,a;y,b) = [5 ab d x ~gf abc A c (x)]6(x~y) is the covari- 
ant derivative in the adjoint representation. It reflects the 
curvature [l2j of the Coulomb gauge field space. In [l2j 
it was proposed to include this curvature in the definition 
of the vacuum wave functional and use, 



$( Wa ,A) = r(yi)$„K,4) 



(4) 



with a = \, which is the usual ansatz for the "radial" 
wave function in curvilinear coordinates. In Ref. fl3| it 
was shown that within the one loop approximation the 
resulting DS equations do not depend on a. We will re- 
turn to this point below. The specific choice a = 1/2 
which has the advantage of eliminating J from the in- 
tegration measure when computing vacuum expectation 
values (vev). 

The DS equations involve expectation values of field oper- 
ators. The gluon two-point function ((Jj'(fe) = 1 — kk/k 2 ) 



dxe lkx (0\A a {x)A b (0)\0) 



5 ab S T (k) 
2u{k) 



(5) 



defines the gluon gap function u>(k), which also relates to 
the single gluon energy. The instantaneous ghost propa- 
gator d(k) (better say ghost two-point correlation func- 
tion) is defined by 



S ab d{k) 



dxe lkx (0\^(x,a;0,b)\0). (6) 



The square of the FP operator —l/(Dd) enters the 
Coulomb potential. The Coulomb form factor f(k) mea- 
sures the ratio of the vev of its square to the square of 
vevs, 



S ab f(k)d 2 (k) = / dxe lkx (0 



Dd 



(x,a;0,b)\0). (7) 



Finally we also consider the expectation value of the cur- 
vature [IH defined by 



S ab d T (k) x (k) 



(8) 

The DS equations for these correlations functions were 
independently studied in (TTJ, [ll[ and [l^, [H| under dif- 
ferent approximation schemes. Below we summarize the 



derivation of these equations and comment on differences 
in their solutions obtained so far. More detailed analysis 
of the equations is presented in the sections that follow. 
The gap equation for the gluon propagator is shown in 
Fig. [T] and follows from minimizing the energy with re- 
spect to the uj a in Eq. (H]), 

ou! a (k) 

In Fig. [TJthc first two diagrams represent the contribu- 
tion from the kinetic energy terms which include trans- 
verse gluon self-interactions from the chromo-magnetic 
energy proportional to B 2 and from chromo-electric en- 
ergy which depends on the curvature and is proportional 
to J~ X EJE. 
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FIG. 2: The DS equation for the ghost propagator, d(k) 
represented by the oval. The gluon line is represented by the 
1/oj propagator. 
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FIG. 1: The gap equation. The diagrams represent vev of 
the Coulomb gauge Hamiltonian. Thin solid line represent the 
derivative with respect to u). The blob represents expectation 
value of one-body operators, and the dashed line represents 
the expectation value of the Coulomb potential. 

The exchange term represents the contribution from the 
Coulomb potential which is determined by the expecta- 
tion value of the Coulomb form factor / and the ghost 
propagator d. Given uj that solves the gap equation © , 
the ghost propagator is computed using Eq. © by ex- 
panding the inverse of the FP operator in powers of the 
gluon field and summing all rainbow-ladder diagrams. 
These are the dominant contributions in the large Nq 
limit. Corrections to the vertex, coupling the transverse 
gluon and two ghosts that appears in the expansion of 
the inverse FP operator were estimated in [l9[ and 
found to be small, of the order of 0(10%). We post- 
pone discussion of rcnormalization to Section IIIII but it 
is worth noting at this point that in Coulomb gauge this 
vertex is UV finite. The diagrammatic representation of 
the DS for the ghost propagator is shown in Fig. [21 
The Coulomb form factor is treated in the same way and 
the resulting DS equation is shown in Fig. [3J Finally, 
to the same one-loop order the curvature is given by the 
loop shown in Fig. U that is determined by the product 
of two ghost propagators arising from the expansion of 
5 2 In J/SASA. 



FIG. 3: The DS equation for the Coulomb form factor f(k) 
represented by the solid circle. The ghost propagator repre- 
sented by the solid oval and the gluon line is represented by 
the 1/ui propagator. 

The resulting DS equations have a general structure, 

F[k, fi(k)] = J dqKlf^k - q), k, q] (10) 

with fi being one of the four two-point function (co,d,f, 
x). While x is merely a definition, u>, d and / obey three 
coupled DS equations. The main interest is in estab- 
lishing if these three equations have a common solution 
and if so what is the momentum dependence of the cor- 
relators, especially in the infrared limit. The bare DS 
equations contain UV divergences and have to be renor- 
malized. Since the variational method does not take into 
account all diagrams contributing at a given order in 
the loop expansion removal of ultraviolet divergences re- 
quires more than just dimensional transmutation of the 
coupling constant. Renormalization of the coupling is 
sufficient to render the ghost propagator finite and thus 
identifies the ghost propagator with the running cou- 
pling. Renormalization of the Coulomb form factor can 
be accomplished by resealing the Coulomb potential by 
a renormalization scale dependent constant. The ghost 
loop is linearly divergent and it turns out that subtrac- 
tion with a finite residual scale dependence is necessary 
to render the gap equation finite. Thus rcnormalization 
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FIG. 4: The one-loop contribution to the FP curvature, x(k). 



of the curvature introduces a scale in addition to the one 
defined by the renormalized coupling. As discussed ear- 
lier additive change of the curvature is equivalent to a 
change in u> a . In particular the inverse gluon propagator 
is given by [H^ 

w(ft)=w a (fc) + (l-2a)x (11) 

For our preferred choice a = 1/2 the wave functional pa- 
rameter coi/2 = ui(k) becomes identical with the inverse 
of the gluon propagator. It follows from Eq. Ijlip that a 
scale-dependent counter-term in x can be absorbed into 
u> a and thus can be interpreted as part of the definition 
of the variational wave functional ansatz. We discuss 
this in more detail in the following section. Finally there 
are quadratic and linear divergences in the gap equations 
which can be removed by the dimension-2 and dimension- 
3 counter-terms proportional to A 2 and {A, II}, respec- 
tively. To summarize, the variational approach intro- 
duced three parameters. One is the scale that can be 
associated with tu a as part of the wave functional defi- 
nition or alternatively it can be associated with \ where 
it represents finite renormalization. The second is the 
value of the coupling constant at some arbitrary value 
of momentum d{k = /i<j), and third is the value of the 
Coulomb form factor at an arbitrary value of momentum 

f{k = n f ). 

The first study of the DS equations in this framework was 
carried out in Ref. [ic| under the approximation X = 0, 
a = 0, so that uj(k) = Wo(fc). The remaining three equa- 
tions were solved analytically under angular approxima- 
tion, 

|fe - q\ -> k0(k - q) + q8{q - k) (12) 

for the integral on the r.h.s. of Eq. (|10[) . Furthermore, 
the solution of the gap equation was approximated by 
w{k) = m g 6(m g — k) + k9(k — m g ). The effective gluon 
mass m g plays the role of the scale parameter discussed 
above. The solutions for d(k) and f(k) were found to 
exist for the renormalized coupling g(fi) = d(k = fx) 
less than a critical value. If the renormalization scale 
is chosen to be m g then g(m g ) < Aiiyj9/{A§Nc) ■ For 



sub-critical couplings both f(k) and d(k) arc finite in 
the limit k — > 0. At the critical coupling d(k) and f(k) 
diverge as 1/ Wk/rrig. Similar behavior of the numer- 
ical solutions without angular approximations was also 
found. For the numerical solutions, however, it could not 
be verified if in the limit g — > g c d(k) and f{k) were 
also IR divergent. This was assumed to be given and ac- 
cordingly an analytical approximation was proposed. In 
[Tl| the effect of the curvature was partially included and 
a numerical sub-critical solution to all three equations 
was found (albeit under some approximations to the gap 
equations). No attempt was made to establish the IR 
limit nor the value of the critical coupling. At this point 
it is worth mentioning that the existence of the critical 
coupling may be an artifact of the variational approx- 
imations and the rainbow-ladder approximation. Func- 
tional integration over the Coulomb gauge fields ought to 
be restricted to the fundamental modular region where 
Det(-Dd) is positive. In our approximation, however, 
in each diagram gauge field integration is unrestricted. 
Thus it is indeed expected that the sum over an infinite 
set of diagrams, (e.g. of the rainbow-ladder series) con- 
verges only in a restricted range of the coupling constant. 

The IR behavior of the Coulomb gauge correlators be- 
yond angular approximation and with full inclusion of 
the curvature was for the first time extensively studied in 
[T3 |. There it was found that in the IR the gap equation 
forces the gluon propagator to have the same momentum 
behavior as the curvature, uj(k) ~ x(fc). Furthermore, 
under the approximation f(k) = 1 it was shown that the 
gap equation and the ghost DS equation admit a solution 
with uj{k) — > oo as k — > 0. Such a solution was not found 
previously in floL [ll[ . In summary, the full set of three 
DS equations have not been solved so far. The differ- 
ences in the various solutions currently available may be 
related to the f(k) = 1 approximation (in the IR limit) 
in pjj or the approximation X = used in [lj|, or other 
approximations made in [ll| . 

In the following we analyze the full set of equations, ex- 
plain why in [12j it was possible to find an IR vanishing 
gluon propagator while it was not the case for the solu- 
tion found in [lol [Tl| , and present a full set of numerical 
solutions to all three equations. 



III. RENORMALIZATION IN THE 
HAMILTONIAN APPROACH 

Each one of the Dyson-Schwinger equations, for the ghost 
propagator, d, the Coulomb form factor /, the gap equa- 
tion for u> and the curvature x requires subtraction of 
the UV divergences. The gap equation and the Faddeev- 
Popov determinant contain power divergences and we 
concentrate on those first. 
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Renormalization of the Faddeev-Popov 
determinant 



Within the approximation we are working the Faddeev- 
Popov determinant 



J (A) = Det(-Dd) = exp (Tr ln(-Dd)\ (13) 
can be replaced by [l3| 



J(A) = exp(- J A X A) 



(14) 



i.e. from the Faddeev-Popov determinant mainly the cur- 
vature x[^](fc) enters. Note that the curvature depends 
on the chosen wave functional and the representation 
Eq. (fT4)) can be used only inside the vacuum expectation 
value (in the state \&) and to the considered order (two 
loops in the energy). From the representation in Eq. (fT4"|) 
it is clear that the counter terms required to renormal- 
ize the Faddeev-Popov determinant, or more precisely its 
logarithm, has to be of the form 



AA 



(15) 



Since the log of the Faddeev-Popov determinant can be 
considered as part of the "action" , we renormalize the 
Faddeev-Popov determinant as 



J -» JAJ = cxp[Trln(-Dd) 

-C X (A) [ dx(A a (x)f 



or by using the representation Eq. (|14|) . we obtain 



JAJ = exp 



- / A( X -C X (A))A 



(16) 



(17) 



Obviously the counter term C X (A) has to be chosen to 
eliminate the ultraviolet divergent part of the curvature 
x(k). Thus the renormalization condition reads in mo- 
mentum space 



X (k)-C X (A) = finite 



(18) 



As usual there is some freedom in choosing the finite 
constants of the right hand side of Eq. (p~5]> . In principle, 
we could just eliminate the ultraviolet divergent part of 
x{k). Since x(fc) has dimension of momentum and fur- 
thermore x(k) is linearly divergent in the ultraviolet, it is 
clear that the ultraviolet divergent part of x(k) is given 
by A • const.. So in principle, we could restrict ourselves 
to remove just the UV-divcrgcnt part of x(fc) by appro- 
priately choosing the counter term C X (A). However, it is 
more convenient to choose C X (A) to be the curvature at 
some renormalization scale /i, resulting in the renormal- 
ization condition 



and in the finite renormalized curvature 
X~(k) = x(k) - x(/«) ■ 



(20) 



It is easy to check that this quantity is indeed ultraviolet 
finite and obviously it satisfies the condition 



X(k = /i) = 



(21) 



By adopting the renormalization condition Eq. (|19p the 
renormalized quantity x(k) Eq. (|20[) depends on the so 
far arbitrary scale fi. By choosing the renormalization 
condition Eq. (|19[) this renormalization scale becomes a 
parameter of our "model" , since it defines the infrared 
content of the curvature x(k), which we keep in the renor- 
malization process, see Fig. [5] 




FIG. 5: The scalar curvature. Note the zero of X (k) a * the 
renormalization scale fj, = 1. 

If we choose for example fx = 0, we chop off the whole 
infrared divergent part of the curvature. 
In Ref. [HI we have chosen a different renormalization 
condition keeping from the ultraviolet divergent quantity 
x(a0 the finite part x'(a0- The renormalized curvature 
then reads 



X (k) = x(k) + x'(M) ■ 



(22) 



This implies that we subtract only the divergent part of 
the curvature keeping fully its finite part. Then we have 
an extra parameter x'(m) °f the theory and the renor- 
malization scale (i is not related to the zero of the renor- 
malized quantity x(k) Eq. (|22|) . Both points of view can 
be adopted: 



1 . Using the renormalization condition of Eq. Q19p re- 
sulting in Eq. (|2ip and the renormalization scale (i 
becomes a physical parameter of the theory. 

2. Removing only the UV-divergent part of the cur- 
vature x(m) keeping its finite part x'(m); which is 
equivalent to choosing the renormalization condi- 
tion 



C X (A) = x(m) 



(19) 



C x (A)=x(//)-x'(M) 



(23) 
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we obtain an extra finite renormalization parame- 
ter x'ifJ-) an d the renormalization scale \i remains 
arbitrary. 



B. Counter terms to the Hamiltonian 

According to the general principles of renormalization 
we should introduce counter terms to the Hamiltonian, 
calculate the expectation value of the energy and then 
choose the coefficients of the counter terms such that all 
divergences disappear from the gap equation. For the 
Yang-Mills Hamiltonian in Coulomb gauge and the wave 
functional at hand it turns out the counter terms required 
are of the form 

AH = Co (A) J dx(A a (x)) 2 

+Ci(A) J dxA a {x) -U a {x) . (24) 

Here the coefficients Cj(A),« = 0, 1 depend on the mo- 
mentum cutoff A and have to be adjusted so that the 
UV-singularities in the gap equation disappear. Note 
also that these coefficients multiply ultralocal operators 
which are singular in quantum field theory 
We will assume the wave functional of Eq. ([4]) with a = 
1/2, as in [l^ . For this wave functional the expectation 
value of the counter term Eq. (|2"4")l is given by 



AE = C (A)-tuS aa J d 3 xu>- 1 (x,x) 

+ id(A) / cPx(A a (x)-iT(x)) u , (25) 



where (. . . ) UJ denotes the expectation value in the Gaus- 
sian wave functional with kernel to — <jJ\/2 and 

ri a (a:) = J^n a (x)J- }2 = n a (x) - - {U a (x) In J) .(26) 
In Ref. uM it was shown 



(A a (x) ■ n°(x)) w = i{A a {x) ■ Q a (x)) u 
i / dx\w~ l (x, x\) {uj(xi, x) — x( x i, x)) 5" 



(27) 



N%-1. In mo- 



where x(k) is the curvature and 5° 
mentum space we find then for the expectation value of 
the counter terms (|251) 



AE = (Nl - 1) V 



Co (A) 



dk 1 



-Ci 



(2tt) 3 w(k) 
dk bj(k) — x(k) 



(2tt) 3 Lu(k) 



(28) 



Taking the variation of this expression with respect to 
w{k) we obtain 



5AE 
5u(k) 



(N 2 C 1) 
-Co(A)- 



V 



(2tt) 3 
1 

J{k~) 



r (\) x(fc) 



(29) 



Note that depends only on the ghost propagator but 
not on the gluon gap function lu, at least as long as x(/c) 
is not yet the self-consistent solution. 
Variation of the expectation value of the energy (without 
counter terms) yields 



6E 
Suj(k) 



N 2 C 



S 3 (0) 



Uk)} , 
(30) 



where 5 3 (0) = j^yi and all other quantities are defined 
in Ref. [l2j]. Adding the counter terms from Eq. (|28[) to 
the energy we obtain the gap equation 



uj 2 (k) 



X 2 (k) = 

k 2 + I° u + I u {k) - 2C (A) - 2d(A) x (fc). 



(31) 



C. Renormalization of the gap equation 



We now turn to the renormalization of the gap equation 
(f3~Tj) . which includes already the counter terms. After the 
renormalization of the Faddeev-Popov determinant, the 
curvature x(k) can b e replaced by the renormalized one 
x(k), Eq. pOj), Note that the Coulomb integral I u (k) 
does not depend on any constant part of the curvature, 
so that x(fc) could have been replaced right away by the 
finite quantity x(^)j Eq. (|2"D|) . Rep lacing x(fc) by x(^) 
and using the relation, sec Ref. 



rtcpl 

m 



uk)=eHk)+2x(k)i^(k), 



(32) 



the gap equation in Eq. (|3ip becomes 

^ 2 (fc)-x 2 (fc) = 

- k 2 + /° + e\k) 2C„(A) + 2 x (k) (/W - d(A)) 

(33) 

The integrals I^J 1-1 ' 2 ^ (k) are linearly (n = 1) and 
quadratically (n = 2) UV-divergent and are defined in 
[12 |. Furthermore, the integral 1° is also quadratically 
divergent but independent of the external momentum. 
We can therefore eliminate all UV-divergences by choos- 
ing Co (A) ~ A 2 and Ci(A) ~ A. In principle, the coef- 
ficients of the counter terms, C(A),i = 0, 1, have to be 
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chosen to eliminate the UV-divergent parts of the quan- 
tities appearing in the gap equation. This means that we 
should choose the infinite coefficients Cj(A), i = 0.1 as 



^+^W) uv _ divergentpart -2C (A) =0 



(34) 



UV-divergent part 



Ci(A) = 0. 



(35) 



Note that the UV-divergent parts of Z^,™ - (fc) are by 
dimensional arguments independent of the external mo- 
mentum k. Technically it is more convenient to choose 
the following alternative renormalization conditions. As 
usual we have the freedom in choosing the renormaliza- 
tion conditions up to finite constants. Given the fact 
that I® is independent of the external momentum and 
the differences 



A/W(ft,i/)=j(")(fc) -/(")(!/) 



(36) 



are UV-finite, we can eliminate all UV-divergences by 
choosing the renormalization conditions 

/°+4 2) (fc = */)-2C (A) = (37) 
l£\k = v)-C{k) = 0, (38) 

where v is an arbitrary renormalization scale, which 
could be chosen to be the same scale \i of the renormal- 
ization of the curvature, but given the fact that with the 
renormalization prcsciption Eq. (|19[) . the scale [i becomes 
a physical parameter, the two renormalization parame- 
ters v and [i need not necessarily be the same. With the 
renormalization conditions Eqs. (|37|38|) the renormalizcd 
(finite!) gap equation reads 

u 2 (k)- X 2 (k) = k 2 + &l£\k lV ) + 2x(fc)A/( 1 >(fc, J ,)(39) 

Assuming that the integrals /f,™ -1,2 '' (k) are infrared finite 
and furthermore that the renormalized curvature x(fc) is 
infrared divergent, the infrared limit of the renormalized 
gap equation is given by 



lim (u>(k) 



x(fc))=A7W(fc = )2 /) 



(40) 



To get the perimeter law of the 't Hooft loop requires 
that the quantity on the right hand side of Eq. (|4"0"1) van- 
ishes. This compells us to choose the, in principle, free 
renormalization parameter v as v = [18l |. With this 
choice the renormalized gap equation becomes 

u?(k) - x 2 (k) = k 2 + A4 2 >(fc,0) + 2 x (fc)A/W(fc,0). 

(41) 

The coupled DS equations then contain the following so 
far undetermined parameters: the renormalization scale 
fi of the renormalization of the curvature and the renor- 
malization constants of the ghost and Coulomb form fac- 
tors, the former one is fixed by the horizon condition. 



D. Renormalization of the ghost and Coulomb 
form factor 



The Dyson-Schwinger equation for the ghost propagator, 
d is given by [l~2| 



1 



d{k) 



1 

9 



Nc 
2 



dl 



1 - (I ■ k) 2 d(l - k) 
UQ) (k - I) 2 



(42) 



and is renormalized by noticing that it corresponds to a 
running coupling. Regularizing the I integral in the UV 
by a cutoff A and replacing g — ► g(A) leads to a finite 
equation for d. It is convenient to renormalize d at an IR 
rather than UV scale which can be done by subtraction, 



Nc f l-{l-kfd(l-k) 

2 J w(0 (k-iy 

(k - M )] . (43) 



d{k) d{fi) 



The Dyson equation for the Coulomb form factor is given 
by 

f ,n 7 A C f u l-{l-k) 2 d 2 {l~k)f{l-k) 

(44) 

Here Zf = Zf(A) is a Coulomb kernel renormalization 
constant and the equation renormalized at an IR scale /i 
is given by 



/(*) = f(p) 



Nc [ dl l ~$- k 



{k-l) 2 



-(fc-p)] . 



(45) 



IV. INFRARED ANALYSIS AT CRITICAL 
COUPLING 

Here, it is shown that, within the approximations made, 
there exists no simultaneous solution to the integral equa- 
tions for d, u), and / at critical coupling for which the 
ghost form factor becomes infrared divergent. We first 
present a simultaneous solution of the integral equations 
for the ghost form factor d(p) and the Coulomb form fac- 
tor f(p) as defined above in Eqs.gS]), (gSJ). The result 
obtained is not compatible with the solution obtained by 
simultaneous treatment of ghost and gluon DS equations 
found in [l9[. We want to check if there are solutions 
of the Dyson equations for u>(p), d(p) and f(p) with the 
ghost form factor d{p) being IR enhanced. To this end, 
we set for p — > oo 

■-HP) = 7-^-,d(p) = -i-J( P ) = ,(46) 



1 

— LO 

2 



(P 2 ) 



(P 2 ) 



(p 2 ) a f 



and k > from the horizon condition. By means of the 
ghost DS equation given in Eq. (|43|) one can relate the 
infrared exponent of the gluon to that of the ghost, 



a z + 2k 



(47) 



8 



in d + 1 dimensional Coulomb gauge YM theory. Thus, 
we eliminate az in favor of k. Using the infrared integral 
approximation, we derive from Eq. (|43|) that 



(p 2 T = (p 2 VAB 2 N c I g (k) 
becomes exact for p — » 0, where 

4"(d-i) r(f-«)r(-«)r(i + «) 

G ( 47r )d/2+l/2 r( d _ 2 At) T(l + f + At) ' 

Furthermore, we can infer from Eq. (|48[) that 
AB 2 NJ g {k) = 1 



(48) 



(49) 



(50) 



has to hold. Simultaneously solving the gluon DS equa- 
tion with u) dominated by the curvature x{p) m the in- 



frared [l^l we find for d = 3 exactly two solutions that 
comply with the horizon condition [l^ . |20| , 



« G <^ 0.398,- 



1 



(51) 



Both of these solutions have been found numerically in 
refs. [l2j and [HI, respectively. 

We now aim at a solution for the Coulomb form factor 
DS equation given by Eq. I|45p. For p — ► 0, one finds 



(p 2 



(p 2 y a fAB 2 N c I f { K ,a f ) 



where 



(52) 



J 



(d-l)« 



T (a/) r (2k) r (d/2 - 2k - 2a/) 



( 47r )d/2 r (d/2 - 2k) T (d/2 

I — 



*/ + l)r(a/ + 2K+l) 



(53) 



From Eq. (|52[) we can see that 

AB 2 NJ f (K,a f ) = 1 



(54) 



has to be fulfilled. The above relation from the Coulomb 
form factor DS equation can now be plugged into the 



J 



relation Eq. (|50[) from the ghost DS equation to find 

Ig{k) = If(K,a F ). (55) 

We now specify the spatial dimension d. For d = 3, Eq. 
([55)1 yields 



1 = 



(-1 + 2 At) cos(tt (a / + 2At))r(4-2a / )r(-2K)r(l + 2a / +4K) sin(7ra/) 
7T (-1 + a/) (3 + 2 k) (-1 + 2 a/ + 4 k) F(2 + 2k) ' 

I 



(56) 



The numerical solution to this equation gives k as a func- 
tion of a/ , as shown in Fig. [5] One can see immediately 
that for any value of a/ , the ghost exponent k yields 



k < 



(57) 



Therefore, recalling the result Eq. (j5"T|) . there exists no 
value of a/ for which all three DS equations are satisfied. 
It is instructive to focus on the case where 



a/ 



2k = 1 



(58) 



since this leads directly to a linearly rising potential for 
static quarks. Plugging this constraint into Eq. ([56]) . we 
find 



-2 k (3 + 2 At) cos(2 7T At) r(-l - 4 At) T(l + 2 At) 
(-1 + 2k) T(-l - 2k) 



(59) 



which has the numerical solution k = 0.245227. Surpris- 
ingly, this result is exactly agreed upon by lattice calcu- 
lations [221 |. where k = 0.245(5) was found. However, one 
should notice that the lattice calculations carried out so 
far in Coulomb gauge and also in Landau gauge use too 
small lattice to give reliable results in the infrared. 
In 2+1 dimensions, the calculations arc equivalent. From 
a simultaneous treatment of the ghost DS equation and 
the gluon DS equation we find for the assumptions that 
the curvature dominates the infrared a unique solution 
for the ghost exponent [2(| , 



(60) 



Note that with angular approximation the result is k = 
1/4. In a recent publication (23[, this value for k was 
confirmed numerically. On the other hand, if we con- 
sider only the ghost DS equation and the equation for the 
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FIG. 6: The ghost exponent function of a/ is shown 

by a solid line for 3 + 1 dimensions. As a dashed line, it is 
indicated for which values a linear potential will arise. 



Coulomb form factor, see Eqs. (|43|) and (|45|) for d = 2, 
one gets by enforcement of the condition Eq. ([55| : 



T(a f ) r(l -a f -2 k) T{k) T(2 + k) 
r(2 - a/) T(-k) 2 r(l+a/ + 2/t) 



1 



(61) 



The numerical solution is shown in Fig. [7J If we require 
the potential to be linearly rising, the solution has to 
obey af + 2k — i . Eq. (ffjTj) then leads to the numerical 
value of k = 0.138. Lattice results do not agree with this 
value t 24|, although they do for 3 + 1 dimensions. 



0.15 




0.05 



FIG. 7: Same as Fig. [6] here in 2 + 1 dimensions. 

It is interesting to note that, leaving the DSE for u aside, 
a restriction on the infared exponents (|5T[) arises from 
self-consistency of the equations for d and / only. This 
restriction also concerns the infrared behaviour of u> via 
the sum rule (|TF)) . Consider a fixed wave functional with 
a kernel ui that yields k > i by its infrared behaviour. 
In principle, it should be possible to calculate d and / 
for such a wave functional. However, due to the restric- 
tion §57} , no solution can be found. We therefore con- 



clude that the Eq. (|57|) must be due to the approxima- 
tions made in the DSEs for d and /. Nevertheless, let us 
stick to the approximation scheme here and investigate 
whether the relaxation of the horizon condition allows for 
a self-consistent solution of all DSEs. 



V. FULL SUBCRITICAL SOLUTIONS 

For the IR analysis given above it follows that the Dyson- 
Schwingcr equations, emerging from the minimization of 
the vacuum energy functional do not admit coupled solu- 
tions with all the functions being IR enhanced. Therefore 
in the following we do not implement the horizon condi- 
tion, but rather chose an infrared finite ghost propagator. 
This allows us to include the ghost propagator fully in the 
equation for the Coulomb form factor. 
Technically, the system is solved in the same way as de- 
scribed in [18|, [21|, except for the following differences. 
First, the finite renormalization constant d^ 1 is imple- 
mented. Then, the infrared extrapolations of the form 
factors no longer can be done using power law ansatze, 
since the form factors are all infrared finite. Instead, we 
approximate the form factors using ansatze of the form 



d ex t{k -> 0) 



1 



kP/B + b 



(62) 



where f3, B and b are parameters and extracted using 
least-squares fitting. One sees immediately that this 
ansatz corresponds to the previous ansatz Bk~@, see (pl6"j) . 
if one sets 6 = 0, and that with a finite b we get for the 
infrared limit of this extrapolating ansatz 



lim d cxt (fc — > 0) 



a— o 



1 



(63) 



which means that the fit parameter b corresponds to the 
renormalization parameter d$ . It can serve as a check 
that the resulting fit parameter b actually turns out to 
be (numerically) equal to the (input) renormalization pa- 
rameter dp . 

If not stated otherwise below, the other renormalization 
parameters are chosen as fi = 1, / M = 1. Physical units 
are not available yet in a strict sense, since the usual 
method of adjusting the Coulomb string tension to its 
physical value can not be employed here, since the re- 
sulting quark-antiquark potential is not linearly rising. 
This is a direct consequence of the finite renormalization 
constant d^ . 

Ansatze of the same form like (|62p are then employed 
for the infrared extrapolations of the other form factors 
x(fc), f(k), oj(k). We chose to use the same ansatze here 
because the numerically obtained values of the form fac- 
tors can be approximated by this ansatz pretty well. Of 
course, every form factor has its own set of fitting param- 
eters. Also x(^) an d w(fc) are fitted separately (instead 
of using a single fit with the combined data, as it has 
been successful in solving the system with the horizon 
condition) . 



10 



A. The search for the critical d(0) 

First, we searched for the maximal value of d(0) for which 
the system still has a solution. If one could carry this 
value to infinity, one had again the horizon condition. It 
is known from the analytical calculations that the sys- 
tem, however, has no consistent solution, once the hori- 
zon condition is enforced. Thus, we expect the system to 
have a solution only for values of d(0) smaller than some 
critical value. 

Results are presented in Figsl8lTT3l In Fig. we show 
the ghost form factor d(k). One finds that the numer- 
ical solutions actually show the desired behavior in the 
infrared, that the infrared fixpoint is at the value of the 
renormalization constant dp . In Fig. [5J we show the 
Coulomb form factor f(k). One finds that with larger 
and larger values of do, the Coulomb form factor f(k) 
becomes more and more enhanced. This has some limit, 
though, as it is not possible to obtain a solution for arbi- 
trary small values of do" 1 . Furthermore, as the Coulomb 
form factor gets more and more enhanced, it becomes 
negative for large momenta. This behavior is not desire- 
able. This fact does not seem to depend sensitively on 
/(/u). The largest d(0) where both f(k) is positive for all 
k and the solution for d{k) is stable has been found to be 
0.02. 



do 1 



100 



10 



d"'(0)=0.l 

Cf 1 (0)=0.027 
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FIG. 9: The Coulomb form factor f(k). 
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FIG. 10: The gluon energy u>(k) and the modulus of the scalar 
curvature, |x(/c)|. The scalar curvature x(k) i s positive in the 
infrared, negative in the ultraviolet, such that in the log-log- 
plot, it seems like it has a pole at its zero. 
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FIG. 8: The ghost form factor d{k). 

In Fig.riniwe show the results for ui(k) and We find 

that these are infrared finite (perhaps with logarithmical 
corrections) , and not diverging like a power law, as it was 
found when using the horizon condition. 
We have stressed so far the new infrared behavior of the 
solutions. As shown in Figs. [T^] and [T2J the ultravio- 
let behaviour is the same as in the case with the hori- 
zon condition implemented. (However, these results are 
hard to compare, since the results with the horizon con- 
dition were obtained using different equations, which use 
the bare ghost propagator in the Coulomb equation, see 
above.) 

Furthermore, we take a look at the product d 2 (k)f(k), 
see Fig. Q31 which is related to the static quark potential 



in momentun space V(k), 

d 2 (k)f(k) = V(k)k 2 . (64) 

Thus, d 2 (k)f(k) ~ l/k 2 corresponds to V(k) ~ 1/fc 4 , 
which after Fourier transform yields a linearly rising 
quark potential. What one finds in the numerical re- 
sults is, of course, that the static quark potential satu- 
rates for large distances, which stems from the fact that 
both d(k) and f(k) are infrared finite, but interesting is 
the feature that for a large momentum range, d 2 (k)f(k) 
actually roughly behaves like 1/fc 2 , which results in a lin- 
ear potential for a large range of distances. We may then 
assume the slope of this linear part as Coulomb string 
tension and introduce a physical scale this way. In Fig. 
1151 the result of V{k) in physical units that were intro- 
duced this way is shown. 
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FIG. 11: Infrared part of the plot shown in Fig. 1101 
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FIG. 12: The ghost form factor d(k), computed as described 
in the text, together with the result described previously, 
where the horizon condition is implemented. Note that the 
Coulomb form factor in the solutions with the horizon condi- 
tion is computed with the bare value of the ghost form factor, 
d = 1, and thus is not enhanced in the infrared. 



B. The search for the critical /(tf) 

In a second calculation, f(fi) was increased to see whether 
a solution where f(k) becomes infrared divergent ("crit- 
ical solution") can be found. However, increasing f(^) 
over some specific value produces a system where Lu(k) 
becomes instable in the infrared. This corresponds to the 
idea that a critical behaviour leads to the non-availability 
of a fully consistent solution of the three coupled equa- 
tions. One can then use an additional degree of freedom 
that arises when one does not set the finite constants on 
the right hand sides of the Eqs. ([57]) and to zero, 
but instead to some finite values. We rewrite the gap 
equation as 



u 2 {k)-x 2 (k) = fc 2 +C + A4 2 )(fc,0) 
+ 2x(fc)(e + A/( 1 '(fc,0) 



(65) 



FIG. 13: C.f. Fig. 1121 but here shown are the gluon form 
factor, ui(k), and the modulus of the scalar curvature, x(k). 
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FIG. 14: The product d 2 (k)f(k). 



where £ and £q are finite parameters. Setting these pa- 
rameters to zero results in the same gap equation ()41l) as 
used before. £g can be used to adjust the value of w(/z). 
In order to keep the (for the 't Hooft loop desired) IR be- 
havior of limfc^o^^) — x(^)) — 0; £ must then assume 
a certain value, 



Co 



2X(0) 



(66) 



Using this additional degree of freedom of tuning £q , one 
can effectively increase iL>(fi), and also the infrared limit 
u)(0). This has the effect that the system becomes stable 
again. However, in the search of the critical /(//), this 
does not help, since the infrared limit of f(k) is decreased 
by the increased infrared strength of u>{k). 
Fig. [TBI illustrates this. We see the solution corresponding 
to Eq. (|4Tj) . with / M = 1. Increasing to 2 produces the 
long dashed line, which shows instabilities in the infrared 
behaviour of ui(k). Increasing £q as well, see the short- 
dashed line, we get stable solutions again, but with a 
smaller (non-critical) value of f(k). 
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FIG. 15: The static quark potential in coordinate space, V(r). 



FIG. 17: uj(k) - x(k) log-over-lin plot. 
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FIG. 16: Solutions for the fully coupled system for different 
/m and 
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FIG. 18: uj(k) - x(k), lin-over-log plot. 



In another calculation, we test the point claimed in the 
introduction that the Coulomb equation is responsible 
for the non-availability of solutions of the fully coupled 
system when the horizon condition is included. To do so, 
we calculate some initial, infrared-finite uj{k) as solution 
of the subcritical system as shown in Fig. [101 We then 
keep this uj{k) fixed and solve the system only for d{k) 
with decreasing dg -1 , and with these d(k) we calculate 
the corresponding f(k). We find that for every finite 
do, both d(k) and f(k) have solutions. With decreasing 
d$ , the d{k) approach a certain asymptotical function. 
In the critical case where we perform the calculation with 
the horizon condition, c£q 1 =0, the ghost equation still 
has a solution, and this solution is equal to the limit of 
the solutions with finite, but decreasing d$ , 

lim ^(fc;^ 1 ) = d(k;d^ = 0) , (67) 

with 

d(fc;0)<oo Vfc>0. (68) 



On the other hand, the f{k) become larger and larger 
with decreasing d$ , and an asymptotical function 
/(fc;^ 1 = 0) that is defined similar to (|67[) does not 
exist. In the critical case where the horizon condition 
is implemented, the Coulomb equation has no solution: 
during the iteration process, the values for f(k) grow 
larger and larger, without limit, to diverge eventually, 
representing the fact that the Coulomb equation has no 
solution in this case. Results of this computations are 
shown in Figs. 1201211 

Finally, we calculated the 't Hooft loop using the solu- 
tions for the form factors presented in this paper. It 
has been shown in [l8[ that the large- i?-behavior of 
the 't Hooft loop is determined by the infrared behav- 
ior of u{k) — x(k)- We give results for this quantity 
u{k) — x(fc), and for the 't Hooft loop itself, in Figs. \Tf[ 
[T9l If we choose our renormalization constants such that 
u(k) — x(fc) approaches for k — > 0, then the 't Hooft 
loop goes to a constant for large R, which is in accord 
with the results given in [l8j |. 
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FIG. 19: The renormalized 't Hooft loop exponent S(R)/R. 



1000 



100 



d" 1 (0)=0.01 
d" 1 (0)=0.005 
d" 1 (0)=0.002 
d" 1 (0)=0.001 
d 1 (0)=0 




1e-04 0.001 0.01 0.1 



100 1000 10000 



FIG. 20: The ghost form factor for different values of dg 1 
calculated with a fixed, infrared-finite uj(k). We find that 
this form factor approaches a limit that corresponds to the 
critical solution with implemented horizon condition. 



VI. SUMMARY 



on the other hand, we did find a self-consistent solution 
where all form factors approach finite values in the 
infrared. This subcritical solution arises if the gauge 
coupling is chosen to be smaller than some critical value. 
Approaching this critical value from below, the Coulomb 
form factor / ceases to solve its DSE. Apparantly, the 
approximation of the ghost-gluon vertex has a different 
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FIG. 21: The Coulomb form factors, calculated with a fixed, 
infrared- finite a>(fc), and the ghost form factors as shown in 
Fig. 1201 This form factor does not have a limit. 



effect on / than it has on d. This explains why in 
previous calculations, critical solutions were found only 
if the DSE for / was ignored. 

The subcritical solution with infrared finite form factors 
presented here have some interesting features. First of 
all, the gluon propagates like massive particle. Secondly, 
with the ghost and Coulomb form factors d and / being 
infrared finite, the heavy quark potential can be found to 
be linear just for a finite range of separation, up to 10/m. 
For phenomenological calculations, the use of these form 
factors may give sensible results, and a further investiga- 
tion would certainly be interesting. 



A non-perturbative solution to the Yang-Mills 
Schrodingcr equation in Coulomb gauge was pre- 
sented. As a variational ansatz, a Gaussian type of 
wave functional was used to minimize the energy. The 
ultraviolet divergences in the three Dyson-Schwingcr 
equations thus obtained were removed by introduction 
of appropriate counter terms in the Hamiltonian and 
in the wave functional, ft was shown that within the 
approximations, no self-consitent solution to all Dyson- 
Schwinger equations can be found that complies with 
the horizon condition. Relaxing the horizon condition, 
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